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Abstract. We study the sedimentation of colloidal gels by using a combination of 
light scattering, polarimetry and video imaging. The asymptotic concentration profiles 
(p{z, t oo) exhibit remarkable scaling properties: profiles for gels prepared at different 
initial volume fractions and particle interactions can be superimposed onto a single 
master curve by using suitable reduced variables. We show theoretically that this 
behavior stems from a power law dependence of the compressive elastic modulus vs 
(p, which we directly test experimentally. The sedimentation kinetics comprises an 
initial latency stage, followed by a rapid collapse where the gel height h decreases at 
constant velocity, and a final compaction stage characterized by a stretched exponential 
relaxation of h towards a plateau. Analogies and differences with previous works are 
briefly discussed. 
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1. Introduction 



Colloidal gels are the focus of an intensive research effort for both fundamental and 
practical reasons. On the one hand, they are model systems to understand the interplay 
between percolation, phase separation and dynamical arrest in systems with attractive 
interactions [U El |3]. On the other hand, colloidal gels are ubiquitous, e.g., in the 
food, drug, personal care and cosmetic industries, where they are often used as a means 
to stabilize a complex formulation against macroscopic phase separation. Mechanically, 
colloidal gels are viscoelastic systems with a predominantly solid-like behavior. However, 
they typically yield under a modest stress, often including the gravitational stress 
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exerted by their own weight. It is therefore not surprising that many studies have been 
devoted to the sedimentation behavior of colloidal gels, revealing a wealth of fascinating 
but yet not fully understood phenomena, such as delayed sedimentation [H El El [71 [8] , 
creep [ll[9l[T0l[TTl[T2], and fracture associated with complex flow patterns |ll[9l[T3]. With 
a few notable exceptions [H El IH], most previous work has focussed on the macroscopic 
behavior of the gels, typically by measuring the time evolution of the gel height, h{t). 
This has prevented a thorough test of the various models proposed to rationalize the 
gel behavior, and in particular of the poroelastic model [Ml [15]. 

To overcome these limitations, we have recently reported in Ref. [16] a full 
characterization of the temporal evolution of the concentration profile, (p{z,t), 
sedimentation velocity profile, v{z,t) and microscopic rearrangement dynamics of a 
colloidal gel. In particular, we have shown that the poroelastic model captures 
remarkably well the evolution of h{t), (f{z, t) and v{z, t). Additionally, a single quantity, 
the compressive strain rate B{t), was shown to control both the macroscopic behavior 
and the microscopic dynamics. Here, we test the generality of the findings of Ref. [I6] by 
studying gels with different particle concentration and interparticle attraction strength. 
We focus on the sedimentation kinetics and on the asymptotic concentration profiles 
as the gels attain mechanical equilibrium, showing that the concentration profiles have 
a universal shape. Using the poroelastic model [HI [15], we show that this remarkable 
scaling behavior stems from a particularly simple relation between elastic stress and local 
volume fraction. The rest of the paper is organized as follows: in Sec. [2] we describe 
the system and the optical methods used to investigate it. The main features of the 
poroelastic model are recalled in Sec. [3l while Sees. 14.11 and 14.21 report our findings 
on the asymptotic concentration profiles and the sedimentation kinetics, respectively, 
before the concluding remarks of Sec. [5] 

2. Materials and methods 

2.1. Sample preparation 

The colloidal particles are spheres of a polytetrafiuoroethylene copolymer (MFA). They 
have a crystalline core, leading to peculiar optical properties that allow one to measure 
accurately the local volume fraction, as discussed below. The data presented here 
are obtained using two different batches, Bl and B2, for which the particle radii are 
slightly different: Ri = 82 ± 3 nm and R2 = 92 ± 3 nm, for batches Bl and B2 
respectively. The high density p = 2.14 g/cm^ of MFA, combined with its lack of 
swelling, allows for a precise determination of the initial particle volume fraction by 
density measurements. The particles are suspended in an aqueous solution of O.IM NaCl, 
to screen electrostatic repulsive interactions, and different amounts of urea were added to 
suppress coherent polarized scattering by matching the solvent and particle refractive 
indices. A nonionic surfactant, Triton XlOO, is added to induce attractive depletion 
interactions [T7] , whose range r ^ 3 nm is of the order of the size of the micelles formed 
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Figure 1. (a): schematic top view of the light scattering apparatus. See the text for 
more details, (b): Typical image of the sample. The gel is the bright, speckled column, 
the black part is the supernatant. 



by the surfactant. In the following we denote by 0tx the concentration of surfactant, 
defined by (f)Tx = Vrx/ (K — Vp), where Vtx, Vp and Vs are the volume of the surfactant, 
that of the particles, and the total volume of the suspension, respectively. A detailed 
description of the phase behavior of this system can be found in |2]. In brief, the system 
displays a metastable liquid-liquid coexistence gap; when the amount of added depletant 
is sufficiently large to drive the colloidal suspension within the coexistence region, as in 
our case, arrested spinodal decomposition leads to the formation of a disordered gel. All 
samples described here are deeply in the coexistence region, far from the spinodal line. 
Because the range of the attractive potential is much smaller than the particle size, our 
system is well described by Baxter's adhesive hard sphere model [H]. It is convenient 
to quantify the strength of the attractive interactions via Baxter's stickiness parameter 
tb- Using methods detailed in [2j, we estimate r ~ 0.01 for the gels of batch Bl and 
tb ^ 0.02 for the gels of batch B2 at <I>tx = 0.11. 

The samples are prepared by mixing a particle suspension and a surfactant solution 
with appropriate concentrations; the resulting suspension is loaded in cells with square 
cross section (either 3x3 mm^ or 5 x 5 mm^) and filled to an initial height Hq = 22.7 mm 
[ho = 32 mm) for batch Bl (B2). The samples are shaken at time tw = and then left 
undisturbed for measurements during the sedimentation process. 

2.2. Light scattering measurements 

The apparatus used for light scattering measurements on systems prepared from the 
first batch is shown in Fig. [1] (a). The sample S is immersed in a transparent water bath 
B to minimize temperature gradients and temperature fiuctuations. It is illuminated by 
a laser sheet LS, of thickness ~ 200 fim, height sufficient to cover the whole gel column, 
and in- vacuo wavelength A = 648 nm. The two lenses LI and L2 form a demagnified 
image of the sample onto a charge-coupled detector (CCD) camera. Their focal length 
is 200 mm and 58 mm, respectively. Two crossed polarizers, PI and P2, are used 
to illuminate the sample with light polarized in the vertical direction while detecting 
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only depolarized scattered light, i.e. light linearly polarized in the horizontal direction. 
Note that this apparatus combines features of both imaging and scattering [19], since 
an image of the sample is formed, but using only light scattered around a well-defined 
scattering wave vector q = AirnX'^ sm{6/2) = 18.5 /xm"^, with n = 1.356 the refractive 
index of the solvent and ^ = 90 deg the scattering angle. A typical image recorded by 
the CCD is shown in Fig. [T](b). The gel appears as a speckled, bright column, while the 
supernatant is black since scattering by the solvent is negligible. Each speckle results 
from the interference of the light scattered by a small volume in the sample, of depth 
equal to the laser sheet thickness and lateral size ~ 50 /im. The speckle size is controlled 
to about 2 CCD pixels by adjusting the range of 9 accepted by the detection optics, 
using the diaphragm D placed in the common focal plane of LI and L2. 

By analyzing a time series of CCD images, we measure four different quantities. 
The gel total height h{t) is simply obtained from the height of the bright column. The 
concentration profiles y^{z,t) are obtained by averaging the intensity / over regions of 
interest (ROIs) of height ~ 100 /xm and width equal to the cell width. We emphasize that 
in general the (polarized) scattered intensity is proportional to both and the particle 
structure factor S{q), thus making impossible a determination of the volume fraction 
based on I without a detailed knowledge of the sample structure. By contrast, our 
particles posses a crystalhne core and thus partially depolarize the scattered light. As 
discussed, e.g., in Ref. [20], the depolarized scattered light can be shown to be simply 
proportional to ip, regardless of the structure. We use an image taken immediately 
after homogenizing the sample to determine the proportionality constant between / 
and ip and to correct for any non-uniformity in the incident beam intensity profile. As 
the gel sediments, the speckle pattern is shifted downwards: using cross-correlation 
methods similar to those used in particle imaging velocimetry [21], we are able to 
measure the full sedimentation velocity profiles, v{z,t), with a vertical resolution of 
about 0.5 mm. Finally, we measure the microscopic dynamics as detailed in Ref. |16j . 
probing in particular particle motion in the horizontal plane, on a length scale of the 
order of ~ i?/1.5. As a final remark, we emphasize that our measurements oi ip, of 
the microscopic dynamics and, to some extent, of v all rely on the assumption that the 
CCD images are formed only by singly scattered photons. Fortunately, the relative low 
average refractive index of the MFA particles (n^ ^ 1.356) allow one to carefully match 
Up using water-based solvents, thereby effectively suppressing multiple scattering. 

For systems prepared from the second batch, we use a simplified version of the 
apparatus shown in Fig. [1], optimized for measuring only h{t) and ^{z,t) [2j. It consists 
of a custom-made light scattering setup, operating at a fixed scattering angle = 90°. 
Selection of the incident and detected polarization of the scattered intensity is made by 
means of two Glan-Thomson polarizers with an extinction ratio better than 1 x 10~^. 
The cell is mounted on a DC- motorized micrometric translator allowing cell positioning 
with a resolution of 0.1 /im and an absolute accuracy of about 3 /xm. A He-Ne laser 
beam is mildly focused in the cell to a spot size w = 46 fj,m, corresponding to a depth of 
focus (Rayleigh range) of about 10 mm, fixing the maximum useful optical path in the 
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cell. The whole setup is enclosed into a removable hood allowing to control temperature 
to better than 0.5° C. 



2.3. Visualization of stresses 

Valuable information on the distribution of the gravity- induced stress in the sample 
can be obtained by polarimetry. The sample is illuminated by a collimated beam of 
white light, issued from a LED source. A demagnified image of the gel is formed onto 
a CMOS camera, working in transmission. Two crossed polarizers are placed before 
and after the sample. Under these conditions, for a stress-free sample essentially no 
light is detected by the camera, because the birefringence induced by the particles is 
extremely weak. Because stress induces birefringence in the particles [22], regions of the 
gel where gravity-induced stress accumulates depolarize the incident light and appear 
as bright in the image (see Fig. [3] below for an example). More precisely [23], when the 
transmission axes of the polarizers are oriented at ±45 deg with respect to the vertical 
direction, stresses along the horizontal and vertical direction are visualized. Conversely, 
by turning the polarizers so that their axes are vertical and horizontal, stresses oriented 
at ±45 deg with respect to the vertical direction can be imaged. 



3. The poroelastic model 

We shortly recall here the main features of the poroelastic model, a popular model 
introduced by Biot [14] and widely used in the literature (see e.g. [I5]) to describe the 
deformation of a gel under its own weight. The system is treated as a continuum medium 
that responds elastically to a compressive deformation. The effect of the solvent back 
flow through the gel is accounted for by introducing a viscous friction term. For the 
sake of completeness, we introduce also a solid friction term, due to the adhesion of the 
gel to the container walls, which is usually not discussed in the literature. 

Typically, the gel sedimentation occurs at a very low Reynolds number and the rate 
of change of v is small; thus, inertia terms can be neglected. Using the reference frame 
shown in Fig. [Tjwith z = the cell bottom, Newton's law for a gel slice of thickness dz 
yields 

dp . da a 

Here, p is pressure, Ap the buoyant density, g the acceleration of gravity, a the elastic 
stress, and L a characteristic length. The term on the l.h.s. originates from the viscous 
drag, while the terms on the r.h.s. account for the buoyant weight, the elastic stress 
and the wall friction, respectively. In particular, L represents Jansen's screening length 
[211 125], such that the gravitational stress in a given horizontal plane S exerted by a 
gel slice located Az above S is divided by a factor exp{AzL~^). This screening is due 
to the redirection of part of the gel weight towards the walls. For L to be finite, the gel 
must "push" against the walls as it is compressed; in other words, the gel must have 
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a positive Poisson's ratio u, such that L = 1{1 — z/)(/iz/)~^ is finite j^H [25] , where / is 
side of the square cell section and fi is Coulomb's solid friction coefficient. We will show 
in Sec. 14.11 that for our gels u ^ and hence L — )■ oo. Thus, in the following we shall 
neglect the wall friction term in Eq. ([1]). 

The conservation of colloidal particles is expressed by the continuity equation: 

|.|W^0. (2) 

The particle velocity v is related to the viscous stress gradient by Darcy's law [26] via 

, = -^|!. (3) 

where k is the volume fraction dependent permeability, 7] the solvent's viscosity, 
g = (1 — ^){vs — v) the fiuid fiux in a reference frame co-moving with the gel, and 
Vs the solvent velocity in the laboratory frame. Since the volume of the suspension 
(solvent plus particles) is conserved, {1 — ip)vs = — v^f , which inserted in Eq. (jS]) yields 

. = ^^. (4) 

7] OZ 

By inserting Eq. (jl]) in the continuity equation, Eq. (|2]), and replacing dp/dz by the 
r.h.s. of Eq. ([1]) with no solid friction term, one finds 

d(p d 
dt dz 



'j^niLp) ( ^ da 



(6) 



■ . , • (5) 

rj \ OZ 

The final equation describing the temporal evolution of the concentration profiles in the 
poroelastic model is then 

dip d 

dt dz \_ 7] \ ' if dz 

where we have used the definition of the uniaxial compressional modulus K 

A»=^|. (7) 

In order to solve Eq. (|6]), the two material functions k{lp) and K{ip) must be 
specified. Insight on K{ip) may be gained by examining the gel behavior for t — )■ oo, 
when the gels approaches mechanical equilibrium. In this limit, v vanishes and so does 
the viscous friction. If wall friction is negligible, at any height z the buoyant weight per 
unit area of the gel column above z is balanced by the elastic response of the gel portion 
below z: 

Apgipdz = cr{z), (8) 



implying that ip must satisfy 

. da da dip _ , 

'^P''^ = -d-z=-d^dI- 
Experimentally, Eq. (|8]) can be used to obtain a{z), provided that <p>{z) can be measured. 
Once both a{z) and p>{z) are known, the volume fraction dependence of the elastic stress 
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is simply obtained by plotting directly a vs (p. As we will show in Sec. I4.H for our gels 
0"((y9) is well modeled by a power law: 

a(^) = A^^ . (10) 

Assuming this functional form, Eq. ([9]) can be easily solved: 

Apg a — 1 



1 
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A a 

where the integration constant Zmax is determined by imposing particle conservation, 
finding 



1 



^max - (^max - /^oo)°"' = /^OV^O K _ 1 \Apg ) ' ^^'^^ 

The determination of the permeability ^{(f ) is less straightforward. In practice, 
we test several functional forms proposed for K{ip) and check whether they are able to 
reproduce the temporal evolution of h{t), the total gel height, as we shall discuss it in 
Seciil 

(13) 



4. Results and discussion 

4-1- Equilibrium concentration profiles 

In this section, we focus on the particle volume fraction and interaction strength 
dependence of the concentration profiles in the long time limit, when no sedimentation 
is experimentally detectable and the gel is very close to mechanical equilibrium. 

For gels of batch Bl and B2, concentration profiles where measured at least 240 
and 500 hours after initiahzing the sample, respectively. We first focus on gels of batch 
B2. Figure [2^) shows (poo{z) = (p{z,t — >■ oo) for three gels, prepared at various particle 
concentrations (po but keeping the same particle interactions (0tx = 0.12). The main 
plot shows that the shape of the concentration profiles does not depend on ipQ-. the 
profiles can be collapsed onto a master curve simply by shifting them along the z axis so 
as to make the top part coincide. Thus, the final height of a gel with a given ip^ may be 
simply predicted using the master curve of Fig. [2] a) and imposing mass conservation. 
Note that this surprising scaling implies that gels prepared at different ipo compress in 
the same way and thus must have the same strength, a counterintuitive result. In fact, 
we anticipate that these gels are formed by the debris of a percolating network, which 
is too weak to sustain itself and collapses almost immediately. These debris must have 
similar structure regardless of ipo, so that the network resulting from their accumulation 
at the bottom of the cell has essentially the same mechanical properties. 

Interestingly, we observe that the concentration profiles are not very smooth. 
Measurements of ipoo{z) for the same gels but using the "beam deflection" method 
described in Ref. [27] exhibit a smoother behavior, probably because in the beam 
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Figure 2. Inset of a): asymptotic concentration profiles for gels with fixed 
0TX — 0.12 and various particle concentrations, as indicated in the labels. Main 
plot of a): same data rescaled onto a master curve as described in the text, b): 
asymptotic concentration profiles for gels with fixed ipo = 0.08 and various surfactant 
concentrations, as indicated by the labels. Lines: fits of Eq. (|III) to the data. 

deflection method the concentration profile is smoothed over a length comparable to 
the cell optical path (5 mm). The noise in ip as measured here may be due also to 
an extra-contribution to the intensity of the depolarized scattered light stemming from 
stress-induced birefringence, which we will show below to be spatially non-uniform. In 
spite of these differences, we emphasize that there is an overall good agreement between 
ip measured by the beam deflection and the depolarized scattered intensity methods [27] . 

The concentration profiles of gels prepared at a fixed particle concentration (po = 8% 
and variable interaction strength are shown in Fig. [2]b). As 0tx increases, interparticle 
attractive forces grow and the gel becomes stronger, leading to a smaller compaction 
under the action of gravity. The change in gel strength is reflected by the fact that 
the shape of the asymptotic profiles is modified. Indeed, in this case the simple scaling 
shown in Fig. [2] a) is no more possible. 

As discussed in Sec. [3], the asymptotic concentration profile is dictated by the (p 
dependence of the elastic stress a, which may be obtained directly from the experimental 
ipoo{z), provided that wall friction is negligible. To investigate whether this assumption 
holds for our gels, we image the distribution of stresses in the gel by polarimetry. The 
inset of Fig. [3] shows the gel column observed under crossed polarizers, where highly 
stressed regions are brighter. In the left image, stresses along the vertical and horizontal 
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Figure 3. Inset: visualization of gravity-induced stress in a gel of batch B2 observed 
between crossed polarizers {(p = 0.08, $tx — 0.12). The height of the image 
corresponds to the gel height, the width of the cell is 5 mm. In the image on the 
left (resp., on the right) the polarizers are oriented so as to visualize stress along the 
vertical and horizontal directions (resp., along directions at ±45 deg with respect to 
vertical). The right image was obtained using an exposure time 1000 times larger than 
that for the left image. Main plot: intensity level along the horizontal line shown in the 
images. The intensity levels have been corrected for the difference in exposure time. 



direction are a priori imaged, although one of course expects vertical stress to dominate. 
The brightness increases towards the bottom of the cell, a first indication that most of 
the gravitational stress is transmitted downward, rather than being redirected to the 
walls. In the right image, the crossed polarizers have been rotated, in order to visualize 
stress oriented at ±45 deg with respect to the vertical direction, which would be non- 
negligible if wall friction was significant. Although some bright patches do appear, we 
emphasize that the right image was taken with an exposure time 1000 times larger than 
for the left one. The intensity profiles along an horizontal line, corrected for the exposure 
time, clearly show that the vertical component of the stress dominates over that at ±45 
deg, implying that wall friction can be neglected for our gels. Interestingly, this is at 
variance with observations for other kinds of colloidal gels and may be due to the 
fact that here the percolating network initially formed is then broken, in contrast to the 
gels of Ref. 

From the asymptotic concentration profiles shown in Fig. [2]b), we calculate crlf) as 
described in Sec. |3] (see Eq. (jH])) and plot the result in the inset of Fig. |4]a). The main 
graph shows that data for all $tx can be collapsed onto a single line in a log- log plot, 
thus indicating that the stress grows with as a power law: a = Aip". The exponent is 
essentially the same for all curves, as shown by the good collapse in Fig. HI By contrast, 
the prefactor A depends on $tx: stronger bonds (larger $tx) yield a stiffer gel. One 
expects A to scale with the interparticle bond spring constant, A ~ Uc/r"^, where Uc 
is the interparticle potential at contact. If the Triton micelles were non- interacting, Uc 




Figure 4. Inset of a): elastic stress as a function of Lp, for gels with i^o = 0.08 and 
$TX as indicated by the labels. Main plot: scaling of cr((p)/A onto a single power law. 
The dotted line has slope 4.08. b) Dependence of the scaling factor A on the amount 
of depletant. 



and thus A should scale with $tx [H], whereas Fig. |4]b) rather suggests A ~ exp[$Tx]- 
Thus, this behavior hints at deviations with respect to the ideal depletion interaction 
induced by a diluted gas of micelles. 

In principle, both A and a may be obtained by fitting (^(v^) to a power law. 
However, we find that a more robust procedure consist in fitting the concentration 
profiles ^{z^t cxd) to Eq. ( ITT]) , the expression derived in Sec. |3] for a power law 
dependence of a vs ip. The resulting fits are shown as lines in Fig. [2]b). We find that 
all the concentration profiles for the gels of Fig. |2] b) can be very well fit using the 
same value a = 4.08. We furthermore check that the fitting parameter Zmax of Eq. (fTTj) 
satisfies Eq. ( fT2i) to within 2%. It is worth noting that, although a must eventually 
diverge as (f approaches the volume fraction of random close packing, no hint of such 
divergence is observed here, nor in Ref. [IE], up to if values as high as 0.45 (batch B2) 
or even 0.55 (batch Bl). 

We check that the law a = A(p°' holds for all gels of both batches, with similar 
exponents: a = 3.6 (resp., 4.1) for batch Bl (resp., B2). Given the power law 
dependence of a vs ip, Eq. (fTTj) suggests that all concentration profiles should collapse 
on the same straight line when plotting ip"'~^ vs {z — -Zmax)/^, where i = A/{Apg) is a 
characteristic length scale that compares the gel elasticity to the gravity pull. Figure [5] 
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Figure 5. Inset: asymptotic concentration profiles for all gels of both batches. Main 
plot: the same data collapse onto a straight line when using reduced variables, as 
explained in the text. Data for gels of batch Bl at ipo = 0.08 (resp., $tx = 0.12) are 
labeled by $tx (resp., ipo). For the gel of batch Bl, ipo — 0.12, ^tx = 0.12. The line 
is the poroelastic model. 



shows that this is indeed the case for all the gels we have studied. 
4-2. Sedimentation kinetics 

The temporal evolution of the total height of the gel, h, exhibits a qualitatively similar 
behavior for all the gels of batches Bl and B2. Denoting by tw the time after filling the 
sample cell, an initial regime where no sedimentation occurs is observed for t„ < td- The 
latency time td grows with the strength of interparticle bonds and ipo, as observed in 
previous works [H El El El E]- During this phase, the speckle dynamics is essentially 
frozen, indicating that a system-spanning arrested network has been formed. For 
td ^ < tc, h decreases steeply and the speckles fluctuate rapidly, indicating that 
the network has failed under gravitational stress. Finally, for t^ > tc the falling debris 
of the initial network have deposited at the bottom of the cell, forming a denser gel, as 
indicated by the dramatic slowing down of the speckle dynamics. This denser gel slowly 
compacts under its own weight: we define t = t^ — tc and describe the gel evolution 
in this third regime using the poroelastic model of Sec. El An example of /i(tw) for a 
gel of batch Bl is shown in Fig. El a). Here, network failure occurs almost immediately, 
td ~ 0. Note that h(t^) is almost identical for gels prepared in cells with square sections 
equal to 3 x 3 mm^ and 5x5 mm^, thus confirming that wall effects are negligible. 
Interestingly, during the initial collapse h{t^) decreases linearly in time, as observed 
also for gels of the second batch prepared at various (fo, see Fig. Elb). This behavior is 
in striking contrast with that reported very recently by Bartlett et al. [8j, who find a 
"compressed exponential" decrease, h(t^ — td) ~ exp{ — [(t^ — td)/Tc]^'^}. This difference 
might be due to the different range of the attractive potential, which is small in our 
case {r/R ^ 0.037) and much larger in Ref. [Ej {r/R ^ 0.62). More experiments will be 



Sedimentation of colloidal gels 



12 



1.0 



0.8 
^° 0.6 
0.4 



0.2 

48 96 144 192 240 

t -tAh) 

w d ^ ' 

1 

10-' 

I 

-« -2 

^ 10 

10-^ 

48 96 144 192 240 

^(h) 

Figure 6. a): time evolution of the total height for a gel of batch Bl with ip = 0.12, 
(/)Tx = 0.12. The cell section is 3 x 3 mm^ (5x5 mm^) for the solid circles (open 
triangles), ha = 22.7 mm. b): initial decay of h during the collapse of gels of 
both batches (same symbols as in c)): in all linear behavior is observed, 

c) semilogarithmic plot of the time varying part of h normalized by ho, during the 
compaction regime. The line is a fit of the poroelastic model to the data for the gel of 
batch Bl. 

needed to elucidate this point. Figure |6]c) shows the time-varying part of h during the 
compaction stage, for various ifQ. The decay of h is not a simple exponential, as shown 
by the curvature of the data in a semilogarithmic plot. Indeed, we find that the data 
can be reasonably well fit by a stretched exponential function, h{t) — ~ exp[— (Ft)'^], 
with /3 = 0.68 — 0.87, depending on sample composition (fit not shown). Although the 
stretched exponential fit works well, it has no obvious theoretical interpretation, besides 
the generic remark that the system exhibits a distribution of relaxation times. A better 
justified fit is provided by the prediction of the poroelastic model, Eq. ([6]), which we 
solve numerically for the gel of batch Bl p6] and show as a continuous line in Figure [6] 
c). As discussed in Sec. El the permeability hi{ip) must be provided in order to solve 
Eq. (Q. By testing various functional forms, including those proposed in the past for 
gels with a fractal structure [12j, we find that only a critical-like law similar to that used 
for suspensions of hard spheres is able to reproduce our results: «(</?) = koV2~^(1 — V^)"^, 
with m = 7. The remarkable success of the poroelastic model is demonstrated by the fact 
that, having determined K{(p) and K,{(f) from fits of (f{z,t — )■ oo) and h{t), respectively. 
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Figure 7. a) sedimentation velocity vs height at various times, ranging from t^/ = 30h 
to tw = 167h from top to bottom, for a Bl gel with (p = 0.12, $tx — 0.12. The line 
is a linear fit to w(2:,tw — 3Gh) up to 2 = 7 mm. b): strain measured over a height z, 
i — \v\/z, for the same data as in a), c): global strain \h\/h (dotted line), z— averaged 
strain e (open circles) and egt as defined in the text, for the same gel as in a) and b). 
The solid lines are \h\/h for the three gels of batch B2 shown in Fig.[6l for which $tx 
increases from left to right. 

the solution to Eq. (|6]) reproduces very well the full time evolution of the concentration 
and velocity profiles, as reported in Ref. [TB] . 

One of the most surprising findings of Ref. [16] is the linear variation of the 
sedimentation velocity v with height z in the gel column. Figure [7] a) shows that 
indeed |f | ~ z, at least for the lower part of the gel column. Close to the top, deviations 
from this linear behavior appear, in that v appears to grow more slowly with z. These 
deviations are more pronounced for longer times. We define a compressive strain rate, 
ifit, as the slope of a linear fit to v{z) in the regime where a linear behavior is observed, 
e.g. up to 2; = 7 mm for data at t„ = 30h in Fig. [7] a). We also define the strain 
rate measured at height z as e{z) = v/z and plot its height dependence in Fig. [7] 
b). This quantity is almost independent of z, with the exception of the very bottom 
and top parts of the gel column, thus confirming that the temporal evolution of the 
full velocity profiles is essentially captured by the evolution of a single parameter, e. 
Additionally, we have shown in Ref. [I6l that e governs also the microscopic dynamics, 
since structural rearrangements occur on a time scale such that Tai = Ey, where 
Ej ~ 0.02 is the typical yield strain beyond which irreversible rearrangements occur. 
In view of the key role played by i, it is interesting to test whether a macroscopic 
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measurement of the global strain rate, \h\/h, reflects accurately the microscopic strain 
rate. In Fig. [7]c) we compare for a gel of batch Bl \h\/h (dotted line) to ^fit (crosses) 
and i, the average of e{z) over the gel height (open circles). In order to reduce the 
impact of the measurement noise on the numerical derivative of h{t), we differentiate 
the stretched exponential flt of h(t) discussed in reference to Figure E]c), rather than 
the data themselves. Figure [7] c) shows that a reasonably good agreement is found 
between these three ways of quantifying the compressive strain rate, thus demonstrating 
that deviations from the linear behavior of v{z) are limited. To check the generality 
of our flndings we simulate the time evolution of gels with different value of K{if ) 
and K,{ip). In all cases we have studied, \v\ ~ 2; in the lower part of the gel column 
while, close to the top, v grows slower with z. Besides iat, i{z) and \h\/h are strongly 
correlated. In particular, we found that \h\/h < e{z) = v/z < ifn ~ A\h\/h with 
1 < A < 1.35 depending on the functional form of K{if) and K{(f), so that e may be 
directly estimated from \h\/h. This observation can be particularly useful when dealing 
with turbid samples, where the velocity proflles can not be measured directly. 

Figure [7] c) shows also the macroscopic strain rate for the same three gels of batch 
B2 as in Figs. [2] a) and El In spite of the variation of (fo, the strain rate at the beginning 
of the compaction stage is of the same order of magnitude for all samples. This is 
probably due to the fact that the structure of the gel is similar, as discussed in relation 
to the scaling of the concentration proflles shown in Fig. |2]a). By contrast, the temporal 
evolution of the strain rate is signiflcantly different, \h\/h decaying faster for the gels 
at lower initial ip. This is at odd with the scaling behavior of the concentration proflles 
reported in Fig. [2^, where it was shown that the top part of the gels at flxed $tx have 
the same concentration proflle. This can be qualitatively understood be recalling that, 
although the top of the proflles for samples prepared at different ipo are similar (see 
Fig j4.ip . gels prepared at higher initial concentration reach a higher volume fraction at 
the cell bottom. Therefore, they display in this region a lower permeability and a higher 
stiffness, which slows down the compaction process. 

5. Conclusions 

We have studied the sedimentation kinetics and the concentration proflles of depletion- 
induced colloidal gels, varying both the particle volume fraction and the strength of 
the interparticle interactions. By using optical methods that combine light scattering 
and imaging, we have been able to gain detailed information on the evolution of the 
concentration proflles and the dynamics of the gels. Coupled to stress visualization 
experiments that rule out any signiflcant role of solid friction on the cell walls, this 
has allowed us to measure quantitatively the volume fraction dependence of the elastic 
response of the gels and to test thoroughly the poroelastic model. In discussing the 
ipo and $TX dependence of the parameters issued from the poroelastic model, it is 
essential to keep in mind that the model can only be applied to the last regime of 
the sedimentation process, i.e. the compaction of a denser gel formed by the falling 
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debris of the initial network that fails under gravity. This allows one to rationalize 
some apparently paradoxical results, such as the scaling of the asymptotic concentration 
profiles obtained at fixed $tx and variable ipo- 

Our work and that of other groups using confocal microscopy [6l [8] show that 
detailed microscopic information on the structure and the dynamics of the gels can 
be very valuable in order to better understand the behavior of colloidal gels under 
gravitational stress. More work along these lines will be required to fully understand 
the microscopic origin of the diversity of behaviors observed for different systems. 
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